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Liouville space search algorithm [Bruschweiler, Phys. Rev. Lett. 85, 4815(2000).] 
utilizes mixed initial states of the ensemble, and has been successfully implemented 
earlier in weakly coupled spins, in which a spin can be identified as a qubit. It 
has recently been demonstrated that n-strongly coupled spins can be collectively 
treated as an n-qubit system. Application of algorithms in such systems, requires 
new approaches using transition selective pulses rather than qubit selective pulses. 
This work develops a modified version of Liouville space search algorithm, which 
is applicable for strongly as well as weakly coupled spins. All the steps of the 
algorithm, can be implemented by using transition selective pulses. Experimental 
implementation is carried out on a strongly dipolar coupled four qubit system. 

1. INTRODUCTION 

An important aspect of information processing is to search a specific data stored in 
a database Classical computer operates on binary codes, any number or symbol is 



represented by combination of bits which can have the value either '0' or'l' 



To 



search an object from N objects, a classical computer must carry out 0(N) operations jll. 
A quantum computer is based on the principle of quantum superposition, which offers the 



opportunity to perform parallel computations 



a 



Grover search algorithm which uses 



quantum superposition, finds the marked object in O(vA^) operations 



m. Most of the 



2 



quantum algorithms, including Grover search algorithm, use a pure state (pseudo pure state 
in NMR) as an initial state HQ, Q. 

An alternative paradigm for computing has been suggested by Madi, Bruschweiler, and 
Ernst, which is based on the mixed states of ensemble, which are represented by density 
operators in the Liouville space Q]. Mixed states describe a statistical ensemble of spins, 
rather than individual spins, so that each element of the ensemble performs part of the com- 
putation, in the same way as each processor in a classical parallel computer l|. Bruschweiler 
proposed a novel search algorithm, that operates on mixed states in spin Liouville space, and 
requires log2N queries to search a single object from N unsorted objects [o^. The advantage 
of Liouville space computation is the use of an input state which is linear combination of 
different classical input states. Protopopescu had apphed Liouville space search to the global 
optimization problem, and achieved exponential speed-up 10, [ll|. Liouville space search 
algorithm of Bruschweiler, considers only weakly coupled spins, in which mixed states are 
represented by linear combination of direct product states of spin polarization operators. 



The algorit 
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13 
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im has been successfully implemented on weakly coupled systems by Xiao et al 
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Most of the NMR Quantum information processing (QIP) experiments have utilized sys- 



tems having indirect spin-spin couplings (scalar J couplings) [16|, ll7|. Since these couplings 
are mediated via the covalent bonds, the number of coupled spins and hence the number 
of qubits is limited to a few qubits. Another approach is to use direct dipole-dipole cou- 
plings between the spins which are larger in magnitude. In liquids the dipolar couplings are 
averaged to zero due to rapid isotropic reorientations of the molecules and in rigid solids 
there are too many couplings yielding broad lines 13]. In molecules partially oriented in 
anisotropic media, like liquid crystals, one obtains partially averaged intra molecular dipo- 
lar couplings with only a finite number of dipolar coupled spins, yielding finite number of 
sharp NMR resonances 13]. However in such cases often the dipolar couplings are large or 
comparable to chemical shifts differences between the homonuclear coupled spins, leading to 



3 



spins becoming strongly coupled. The method of spin selective pulses and J evolution used 
in NMR-QIP of weakly coupled spins (liquid state NMR-QIP), leads to complications in 
strongly coupled spins because of U te.n. in tKe Ha,nnton,a„ Howeve. it Ka. been 
demonstrated that the 2" non degenerate energy levels system can be collectively treated 
as an n-qubit system, similar to the case of quadrupolar nuclei oriented in liquid crystal 



19 



20 



2l|. The manipulation of the states of the qubit are then carried out by using tran- 



sition selective pu 
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2Q|, 
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ses. Several gates and algorithms have been implemented in such systems 



24 



25|. 



In this work, we generalize the Liouville space search algorithm such that it can be 
implemented in strongly as well as weakly coupled spins. Experimental implementation is 
carried out on a four qubit system obtained by strongly dipolar coupled protons of 2-chloro, 



iodo benzene oriented in ZLI-1132 liquid crystal [2^, l25|]. In section (2) we explain Liouville 
space search algorithm, in section (3) we explain the modified version, section (4) contains 
the experimental implementation, and conclusions are given in section (5). 



2. LIOUVILLE SPACE SEARCH ALGORITHM 



Liouville space search algorithm finds a marked object in 2" unsorted objects, which 
correspond to 2" eigen states [9|. The eigen states in weakly coupled spins, are same as 
product states obtained from individual spin sates, and the corresponding density matrices 



are obtained from direct product of spin polarization operators 



18| 



= \al3a a) = \a)\j3)\a) \a) 

a = \^P){^\=I'^I^I^ i:, 



where /"/g /f" is the direct product of polarization operators J^, jf , defined as. 



4 




It = \ak){ak\ = -{lk + 2h,)-- 

lt = \(3u)m = \{tk-2h-:)= \ y \. (2) 

21 kz is the Pauli matrix and 1^ is unity operator of the subspace of spin 1^. In Eq. 
([1]) a has the dimension of 2" and it's matrix elements in the eigen basis (also known as 
Zeeman basis or product basis) are all zero except for one diagonal element belonging to the 

state \a[3a a) with coefficient 1, which represents the population of that state (Eq. [1]) 

as one. |a) and are assigned to the computational basis states |0) and |1), thus Eq. ([T]) 
can be written as, 



i) = |oio 0) 

a = \^P){^\ = I'^nI! C (3) 

The algorithm aims to search a single object in = 2" unsorted objects which are 
represented by 2" basis states of 'n' qubits, which we call as work qubits. The density 
matrices of these states are represented in the Liouville space as, 



-"n 1 1 



(Tl = I'J^II J°= |0i0203....0„)(0i0203....0, 

(72 = I^I^ll J°_i/^= |0i0203....0„_iU(0i0203....0„_il„| 

a2n = /;^ = |lil2l3....1„)(lil2l3....1„|. (4) 

Let the marked state is am (m=l,2,....,2" ). The algorithm uses an extra qubit known as 
ancilla qubit Jq, and finds the marked object in 'n' (= log2N) queries. 
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(i) Preparation of initial states: 'n' initial states, (k=l, 2,...,n), are prepared in 
'n' different experiments, pf„ is given by [o], Q , 



Pin — -'o-'fc 



= /0II1I2 tk^Jltk+l In] 

= I>i+^2+ + = I'oal (5) 

In Eq. ([5]), ancilla qubit and A;*'^ work qubit are in state |0), whereas each of the remaining 



lus pf„ has equal and non-zero populations (each of 



[n — 1) qubits are in state |0) or |1). T 

magnitude 1) in 2"~^ energy levels j^, I2I. In other words, all the matrix elements of pf„ 
are zero, except for 2"~^ diagonal elements, each of which represent populations of value 1. 
The initial state /3^can be prepared in weakly coupled spins by using spin selective pulses 



and J-e volutions 



12 



13|. 



In Eq. ([5]) since fc*'' work qubit of each af, is already set to \0k), if crf„ contains the 
marked state a-m, then it implies that A;*'^ work qubit of the marked state is |Ofc). In other 
words (Jm= Ii^^l2^^ ■■■■Ik-i^iilk+i----'^^^^ ■ '^in ^ocs uot coutaiu (Tm, then k^'^ work qubit of 
the marked state must be that is 0"^= /|'^^/2^\.../°^\lJ^/°^\..../n''^ 0, Q- 

(ii) Oracle: Oracle is a black box, contents of which are unknown to the reader. Unitary 
operator (U) of oracle, inverts the states I^am and lQ<7m- The resultant state is then given 

by, 



p^ = f/(ptj) = U^I^crfj = Jocrf^, if does not contain the marked state am, h =^ H 

if contains the marked state am, h =^ Ik- (6) 
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As one can see in Eq. ([6]), 'U' is operated simultaneously on all the qubits of the ensemble, 
(iii) Measurement: The ancilla qubit (Jq) is detected by applying a (vr/2)° pulse on 
[121]. The spectrum of Iq contains 2"~^ peaks corresponding to 2"~^ states of af^. In 
Eq. (EI), if = then all the peaks are of positive intensity, which indicates that the 



12|. If 



k*^ qubit of the marked state is in state |1) or crm= ^i^^ ^2^^ ■■■■Ik~i^l.'^k+i----'^n^ [9 

'^IcTrr, — lolo^in ~ d"]) then the spectrum contains a negative peak corresponding to 
am, whereas remaining (2"^^ — 1) peaks are of positive intensity, which indicates that am= 



12 



.0/1 .0/1 .0/1 yo rO/1 rO/l 
-'l -'2 •■••-'fc--l-'-k-'fc+l----'" 

A sequence of 'n' experiments with input pf„ = IqI'^, k=l, 2, n, are recorded. In the 
k^^ experiment if all the peaks are positive then the marked state contains I\, otherwise it 
9|. Thus in 'n' experiments one can find the states of all the work qubits, Ji, /2, , 



is 4° 



In of the marked state. 



3. GENERALIZED LIOUVILLE SPACE SEARCH ALGORITHM 



Liouvi 



spins 



12 



le space search algorithm has been implemented in NMR by using weakly J-coupled 



13| . Various steps of the algorithm are implemented by using qubit (spin) selective 



pulses and evolving the system under J-couplings 12, llSl] • K is difficult to implement this 
algorithm in strongly coupled spins, since the spins are not individually addressable in the 
spectrum. Also it is difficult to manipulate the time evolution under a desired Hamiltonian. 
In other words eigen states are not the product states of individual spins. However an n- 
coupled system can be collectively treated as an n-qubit system, by treating the 2" eigen 



states as computational basis states 



21 



22 



2J]. In such systems it is difficult to apply 



qubit selective pulses, but one can still manipulate the states by applying transition selective 
pulses. Thus the original algorithm explained in section (2), is modified in such a way that 
various steps of the algorithm can be implemented by using transition selective pulses. The 
modified algorithm is implemented here in strongly coupled spins, but can be implemented 
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in weakly coupled spins as well. 

Consider an (n+l)-qubit system (weakly or strongly coupled), where the 2""'"^ eigen states 

are labeled as |OoOi....O„) |loli....ln), where the subscripts '0' and 1, 2...., n represent 

ancilla qubit and 'n' work qubits respectively. The search space contains 1^ states of 'n' 
work qubits, whose density matrices in 2"- dimensional subspace of 2"^^- dimensional eigen 
basis, are given by. 



Ox = = IO1O2 0„)(0i02 0„| 

<J2 = \^P2){H = IO1O2 I„)(0i02 1„| 

O-gn = |V'2")(V'2"| = II1I2 In)(lll2 ln|- (7) 

In Eq. ([7]), the matrix elements of cTj [i = 1,2, ....,2") are zero, except for one diagonal 
element of value 1, which corresponds to the state ipi- The marked state is represented by, 

Om = |^m)(^m| = \x1X2 Xn){xiX2 x„|,where m=l, 2,..., 2^" and Xi=0 or 1. (8) 

(i) Preparation of initial states: 'n' initial states (k=l, 2, ...,n) are prepared in 'n' 



different experiments using the method of POPS (pair of pseudopure states) [l9|. POPS(i, 
j) contains populations only in states \i) and given Iry. 'p' and '-p' respectively, where 
the value of 'p' depends on the details of preparation [19]. POPS(i, j) can be achieved 
by inverting the populations of a pair of states \i) and \j) and subtracting the resultant 
population distribution from that of equilibrium state 19|. 

In each initial state p^^, ancilla and fc*^ work qubits are in state |0), whereas each of 
the remaining (n-1) qubits are in state or 1. Thus the initial state contains non zero 
populations in 2"~^ states, also for every state of population pf there exist another state 
of population of —pf. In other words, pf^ is a sum of {2^^^/2) POPS, thus it represents a 
highly mixed state, pf^ is prepared by subtracting pk from peq (equilibrium populations). 
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where pk is obtained by inverting the populations of (2" ^/2) pairs of states, and subtracting 
the resultant population distribution from that of equilibrium state, pf„ is given by, 

k 

Pin = Peg " Pk, 

= Pt|0o0i02....0fc_i0k0fc+i....0„)(0o0i02 Ok 0„| 

+P2|0o0i02....0fc„i0k0fe+i....l„)(0o0i02 Ofc_iOkOfc+i 1„| 

+ +p^„-i|Oolll2....1fc-lOklfc+ l....l„)(Oolll2 Ifc-lOklfc+l In. I 

= P?|OoV^l)(OoAl +P2\00^2){00^2\ + +P2'"-i|Oo^2"-i)(OoV^2"-i|, (9) 

where p^,P2----P2n-i represent the populations in the respective levels, and for every pf 
there exists a -pf, thus pf„ represents a sum of (2""^/2) POPS. In Eq. ([9]) the k^'^ work 
qubit is set in state \0k), hence if pf„ contains the marked state \iprn) = \xiX2..--Xk---.Xn), it 
implies that the fc*^ work qubit of lipm) that is Xk = 0, whereas if pf„ does not contain \ipm) 
then Xk = 1- 

(ii) Oracle: The unitary operator of the oracle acting on pf^, flips the states |Oo^/'m) and 
jloV'm)- 'U' can be implemented by applying a selective (vr) pulse on transition \Oo'ipm) 
lloiprn)- The final state p'j is given by, 

p'j = U{p\^) = pt).,if Pin does not contain the marked state \il)m) =^ Xk=l 

if contains the marked state \ipm) =^ Xk=0. (10) 

(ill) Measurement: A multi frequency (7r/2)j, pulse is applied on 2" unconnected 

transitions (ancilla qubit transitions) corresponding to 2" pairs of states \Ooi1i2 in) 

\loi1i2 in), where ii, 12, ■■■■in = or 1 , followed by the detection. In case of weakly cou- 
pled system with conventional labeling (that is the states \a) and \(3) of each spin, represent 
|0) and |1) respectively), these 2" transitions belong to Iq qubit. Hence the read-out pulse 
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is selective {T^/'2)y pulse on Jq qubit. Thus in case of weakly coupled spins, measurement 
strategy is same as that of the original algorithm [o], [section 2]. 

In Eq. f lTUl) . if the final state Pj=/5f„, then the MF pulse converts the populations 

Pi|Oo'0i)(Oo^5'|, P2|Oo'02)(Oo^2l> P2"-i|Oo^2"-i)(Oo^2"-il' ^""^ cohereuccs. Since for 

every p\ there exists -pf, number of positive transitions (associated with positive pf) is 
equal to number of negative transitions (associated with negative p^). If contains the 
marked state (corresponds to second equality of Eq. [10]), then positive and negative pop- 
ulations (other than pJ^lloV'm) (loV'ml) are converted in to positive and negative transitions 
respectively, whereas a positive or negative p^|loV'm)(lo'^m| gives a negative or positive peak 
respectively. Thus, if contains the marked state, then spectrum contains (2"~^/2) ± 1 
positive peaks and (2"^Y2) =F 1 negative peaks, where + and - respectively corresponds to 
negative and positive p^. 

In summary, in the k^^ experiment with initial state p^„, if the spectrum of final state 
contains equal number of positive and negative peaks, then the A;*'* work qubit is in state |1) 
or Xk = 1. Whereas if the spectrum contains unequal number of positive and negative peaks 
then Xk = 0. This procedure is repeated for each of the 'n' initial states pf„ (k=l,2...n), then 
one obtains the marked state a^. = \x1X2 x„)(xiX2 a;„|. 

It is to be noted that, since pf„ = Peq—pk (Eq- E]), implementation of 'U' and measurement 
are performed in different experiments for peq and pk (k = 1, 2,....,n). Now the spectrum 
of pj is obtained by subtracting the final spectrum of pk from that of peq- Thus the total 
number of experiments is (n+l). 

4. EXPERIMENTAL IMPLEMENTATION ON A STRONGLY DIPOLAR 

COUPLED 4-QUBIT SYSTEM 

The system chosen is, four protons of 2-chloro, iodo benzene oriented in ZLI-1132 liquid 
crystal [22]. The four protons are strongly coupled to each other through dipolar couplings. 
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and contains 16 eigen states. The total Hamiltonian of the system contains, chemical shift 
(He), dipolar (Hod), and J-coupling (Hj) terms, given by. 



H = H, + Hdd + Hj = Y^ uJi, + ^ [27rA,(3/i./,. - h.Ij) + 27r Jy/,./,]. (11) 

i=l i,j=l,andi<j 

Experiments are carried out on AV-500 NMR spectrometer. Figure ([T]) shows the 
ID spectrum obtained by applying a hard {Tx/2)y pulse. In order to perform the Liouville 
search, one needs to know the pair of states from which each transition originates or in other 
words, we need to know the complete energy level diagram of the spin system. The energy 
level diagram of this system, has been obtained by Mahesh et al, using the connectivity 
information from Z-COSY experiment The pulse sequence of a Z-COSY experiment 

is, 90- ti - a - T - (3 - 12-, where only the longitudinal magnetization is retained during the 
interval r by either phase cycling or by a gradient pulse 3, 26]. The 90" pulse converts 
the equilibrium z-magnetization into coherences, which are frequency labeled during the 
period ti. The small angle a pulse ensures that each cross-section parallel to uj2 from the 
resulting 2D spectrum is equivalent to a one-dimensional (ID) experiment in which the peak 
corresponding to the diagonal is selectively inverted 26[ . The directly connected transitions 
to the inverted transition are finally measured in the linear regime by a small angle j3 pulse, 
igure ([2]) shows the energy level diagram of this system, obtained from Z-COSY experiment 



22 



where each of the transitions of Fig. ([T]) are assigned to a pair of energy levels. 
The 16 eigen states (Fig. [2]) are labeled as computational basis states of a four qubit 

system, given by, IO0O1O2O3), IO0O1O2I3), II0I1I2I3), where the subscripts '0' and 1, 2, 

3 represent ancilla qubit and three work qubits. It is to be noted that the labeling of eigen 
states shown in Fig. ([2]) , is different from that of ref . [22I] . The labeling in our case follows two 
strategies: (i) Find eight unconnected transitions, and label the eigen states such that these 
transitions belong to ancilla qubit. In other words, each of the eight transitions is associated 
with a pair of states |Oo^l^2^3) -v^ |lo^i^2^3), where ziZ2^3 = 000,001, ,111; the states 



11 



|Oo^i^2^3) and |lo'^i^2^3) respectively represent the lower and upper states of the transition. 
However in this system we could find only seven unconnected transitions, represented by 
dark lines in Fig. ([2]), the transition between the states IO0I1I2I3) II0I1I2I3) is not 
observed due to low intensity, (ii) The initial states pf„ (Eq. [13]) contain the eigen states 
|0000), |0001), |0010), |0100), |0011) and |0101), hence to prepare pf„ each of these states 
should not be isolated from the rest of eigen states, for example the labeling of any of these 
eigen states can not be interchanged with that of |0111) since this state is isolated (Fig. [2]). 
Following these two criteria there exists many possibilities of labeling these eigen states, the 
present being one of them. It may be pointed out that the second strategy is not necessary, 
if none of the states are isolated. 

Liouville space search is now performed, following the procedure given in section (3). The 
search space contains eight states of the work qubits, given by. 



(Ti = 




= |0l0203)(0i0203|. 


(72 = 


= \i^2){'ip2\ = 


|0l02l3)(0l02l3| 


0-3 = 




= |0il203)(0il203|. 


(T4 = 


= \^pA){i^A\ = 


|0il2l3)(0il2l3| 


= 


1^5) {^bl 


= |ll0203)(ll0203|. 


0-6 = 


= \ = 


|ll02l3)(ll02l3| 


ar = 




= |lil203)(lll203|, 


and 


as = \^s){^i 


j| = |lil2l3)(lll2l3| 



(12) 

Let the marked state is ctq = |'^6)(V'6|- 
(i) Preparation of initial states: 

Three initial states pf„ (k=l, 2, 3) are required since there are three work qubits. From 
Eq. ([9]) we can write these initial states as. 
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P]n = p}|OoOl0203)(OoOi0203| +P^|OoOi02l3)(OoOl02l3| 

+p^|OoOil203)(OoOil203| +P^|OoOil2l3)(OoOil2l3| 

pI = P?|OoOl0203)(OoOi0203|+P^|OoOi02l3)(OoOi02l3| 

+P3|Ooll0203)(Ooll0203| +P^|Ooll02l3)(Ooll02l3| 

pL = P?|OoOl0203)(OoOi0203|+P^|OoOi02l3)(OoOi02l3| 

+p^|Ooli0203)(Ooll0203| +p^|Ooll02l3)(Ooll02l3|. (13) 

In each experiment, for every there exists another population which is -pf . 

A density matrix pi is prepared from p^g (Fig. [2]), by applying a selective vr pulse on tran- 
sition (23) which interchanges the populations of IO0O1O2O3) & IO0O1I2O3), and a cascade of vr 
pulses on (30), (9) and (30) which interchanges the populations of IO0O1O2I3) & IO0O1I2I3). 
Now p}^ is obtained by subtracting pi from peq, the resultant populations of p}^ are shown 
in Fig. Q, which can also be written as, 

pL = Peg-pl = |OoOl0203)(OoOi0203| + 2|OoOi02l3)(OoOi02l3| 

-|OoOil203)(OoOil203| - 2|OoOil2l3)(OoOil2l3|. (14) 

Figure (j4^) shows the spectrum of pj„ (Fig. [3]), obtained by subtracting the spectrum of 
pi from Peg, where each of the spectra of pi and peq were recoded by applying a 5° pulse. In 
Fig. ([3]), the population difference corresponding to transitions (29), (24), (23), (27), (30), 
(25) and (9) are positive, while those of (16), (22), (28), (14) and (1) are negative. This is 
confirmed from the spectrum of Fig. (jl^), the transitions (14) and (16) are not visible due 
to low intensity. 

The is given by. 
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pI = peg -P2 = |OoOi0203)(OoOi0203| + IO0O1O2I3) (O0O1O2I3I 

-|Ooll0203)(Ooll0203| - |Ooll02l3)(Ooll02l3|, (15) 

where p2 is obtained from p^q (Fig. [2]), by applying six selective vr pulses on transitions 
(24), (30), (21), (5), (21) and (30) respectively, where the first pulse inverts the populations of 
IO0O1O2O3) & IO0I1O2O3) and next five pulses invert the populations of IO0O1O2I3) & IO0I1O2I3). 
Figure ([5]) shows the population distribution of pf^ and the corresponding spectrum is shown 
in Fig. (jib). The positive and negative transitions in Fig. (jib) are in accordance with the 
population distribution of Fig. (j5]) . 

Similarly one obtains, 

P\n = Peq " PS = | O0O1O2O3) (O0O1O2O3 | + | OqOi I2O3) (OqOi I2O3 | 

-|Ooll0203)(Ooll0203| - |Oolll203)(Oolll203|, (16) 

where p3 prepared from p^q (Fig. [2]), by applying two selective vr pulses on transitions (24) 
and (22) which invert the populations of IO0O1O2O3) & IO0I1O2O3) and IO0O1I2O3) & IO0I1I2O3) 
respectively. Figure (jit) shows the spectrum of p\^ which reflects the population distribution 
of Fig. (jej). 

(ii) Oracle: To search the state \iPq) (Eq. [12]), the unitary operator (U) of oracle, 
interchanges the states jOo^/'e) and jlo^'e)- Here U is achieved by applying a selective vr pulse 
on transition (18). The resultant states pj-, and p^ are given by. 
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p) = u{pl) = u{p,,) - u{p,) = pI 

P) = Uipl) = U{p,,) - Uip2) 

= f/(|OoOi0203)(OoOi0203|) + f/(|OoOi02l3)(OoOi02l3|) - f/(|Ooll0203)(Ooll0203|) 

-f/(|Ooll02l3)(Ooll02l3|) 
= |OoOi0203)(OoOi0203| + |OoOi02l3)(OoOi02l3| - |Ooll0203)(Ooll0203| 

-(|loll02l3)(loll02l3| 

P/ = Uipl) = Uipe,) - U{ps) = pI. (17) 

In Eq. (fT7|) . U converts the state IO0I1O2I3) of in to II0I1O2I3), whereas p}^ and pf„ 
are not effected by U, since neither of these states contain either IO0I1O2I3) or II0I1O2I3). 

(ill) Measurement: Measurement requires a multi frequency (vr/2)y on eight uncon- 
nected ancilla qubit transitions of the final state (pj) followed by the detection. The eight 
transitions of MF pulse correspond to eight pairs energy levels, \Ooi1i2i3) 1 10^1^2^3)1 where 
^i,'^2,'^3 = or 1. The transition between the energy levels IO0I1I2I3) and II0I1I2I3) is not 
observed in the Z-COSY assignment of Fig. ([1]) 2^, due to low intensity. However this 
transition is not required, since it any way gives rise to zero intensity peak, because the 
final states pj after the oracle of search sate {i/jq) (Eq. [T7|) and also for remaining search 
states, contain zero populations in IO0I1I2I3) and |lolil2l3)- Hence the multi frequency 
(7r/2) pulse is applied on remaining seven transitions, (1), (8), (10), (18), (27), (28), (29) 
(represented by dark lines in Fig.s[3], [5], [6]). The duration of MF pulse is 70 ms, which is 
obtained by modulating the Gaussian pulse with seven harmonics corresponding to these 
seven transitions, the phase of each harmonic is 'y'. The amplitude of each harmonic is 
adjusted such that the pulse gives maximum intensity for each transition, in other words, 
the MF pulse corresponds to a (vr/2)y pulse on each of the seven transitions. The spectra 
of pp p'j and pj are given in Fig.s[7| (a), (b) and (c) respectively. In Fig.s ([7^, c) number of 
positive peaks is equal to number of negative peaks, hence the first and third qubits of the 
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marked state, are in state |1). Whereas in Fig. ([Tb), number of positive peaks is not equal 
to number of negative peaks, hence second work qubit is in state |0). Thus the marked is 
IliOsls). 

It may be pointed out that, since the ancilla qubit transition |0111) is not 

observed (Fig. [HE]), one can not perform the Oracle operation (U), to search the state |111) 
with the present labeling scheme. 



5. CONCLUSIONS 



In order to increase the number of qubits, one has to exploit the dipolar couplings among 
the spins, in which case the spins are often strongly coupled to each other. Unlike weakly 
coupled systems, strongly coupled systems can not be directly used for implementing quan- 
tum algorithms. In this work, we generalize the Liouville space search algorithm, such that 
it can be implemented in weakly as well as strongly coupled systems. Experimental imple- 
mentation is carried out on a strongly dipolar coupled four qubit system. All the steps of 
the algorithm are implemented by using transition selective pulses. 
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FIG. 1: H equilibrium spectrum of 2-chloro, iodo benzene oriented in ZLI-1132 liquid crystal, 
recorded at room temperature (300k) on a 500 MHz NMR spectrometer. The four protons of the 
molecule are strongly dipolar coupled to each other. Various transitions are numbered in increasing 
order from left to right (decreasing frequency). 
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FIG. 2: Energy level diagram of oriented 2-cliloro, iodo benzene, obtained by investigating the 
connectivity of various transitions of Fig. [Uby using a Z-COSY experiment [22]. The equilibrium 
populations {peq) of various levels (represented by bold numbers) are proportional to the Zeeman 
energies of the protons, schematically represented by the energy level diagram. The 16 eigen states 
are labeled as basis states of a four qubit system. The strategy adopted for this labeling scheme 
is explained in the text. The transitions represented by dark lines, correspond to ancilla (zeroth) 
qubit. 
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FIG. 3: Energy level diagram of Fig. ([2]), with the populations of various states corresponding to 
pjj^ (Eq. [T^ . The transitions represented by dark lines, correspond to ancilla qubit. 
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FIG. 4: The spectra of (a), (b) and (c) respectively represent the population distributions of pj^ 
(Fig. [3]), (Fig. [5]) and (Fig. [6]). (a), (b) and (c) are respectively obtained by subtracting the 
spectra (recorded with 5° pulse ) of pi, p2 and from equilibrium (peg) spectrum (also recorded 
with a 5° pulse). Preparation of pi, p2 and ps are explained in the text. Some small intensities 
such as 14 & 16 in (a), are not observed/marked. 
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FIG. 5: Energy level diagram of Fig. ([2]), with the populations of various states corresponding to 
p?„ (Eq. [T5]) . and the transitions represented by dark lines, correspond to ancilla qubit 



24 







I 1 1 i> 




FIG. 6: Energy level diagram of Fig. Q, where the populations of various states correspond to 
(Eq. [T6|) . The transitions represented by dark lines, correspond to ancilla qubit 



25 



(b) 



T 




29 


1 


1 

2S 


6000 5000 4000 3000 2000 1000 

2 

27 

J 1, 


1 

Hz 

9 





6000 5000 4000 3000 2000 1000 Hz 



6000 5000 4000 3000 2000 1000 Hz 

FIG. 7: Implementation of Generalized Liouville space search algorithm on a dipolar coupled four 
qubit system. The spectra of U(/9e(j)) U(/9i), U(/92) and U(p3) (Eq. [T7|) are first obtained by 
applying a multi frequency pulse on ancilla qubit transitions (1), (8), (10), (18), (27), (28) and 
(29). Then the spectra (a), (b) and (c) are respectively obtained by subtracting the spectra of 
U(/3i), U(/52) and U(p3) from U(pei3) spectrum, (a), (b) and (c) thus respectively represent the 
ancilla qubit spectra of pj, and (Eq. Since /jj- = and = p?^, the signs of ancilla 

qubit transitions of (a) and (c) are respectively equivalent to that of Fig. s (H^) and ([Ih), whereas 
in (b) of the sign of transition (8) is changed from that of p^^ Fig. (Hb), due to the oracle 
operation. Thus the spectra of p^ and pj (a,c) contains equal number of positive and negative 
peaks, whereas the spectrum of pj (b) contains unequal number of positive and negative peaks. 
Thus the search state is |101). 



